Journal of Soft Computing in Civil Engineering 2-1 (2018) 71-84 


3A 
Cc 
O 
= 


press 


Contents lists available at SCCE 


Journal of Soft Computing in Civil Engineering 


Journal homepage: www.jsoftcivil.com 


Flexural Analysis of Deep Aluminum Beam 


A. Dahake'’®, P. Kapdis’, U. Kalwane’, U. Salunkhe’ 


1. Associate Professor and Head, Civil Engineering Department, Maharashtra Institute of Technology, Aurangabad 


(M. S.), India 


2. Civil Engineering Department, Shreeyash College of Engineering and Technology, Aurangabad, (M. S.), India 
Corresponding author: ajaydahake @ gmail.com 


5) _https://doi.org/10.22115/SCCE.2018.49679 


ARTICLE INFO 


ABSTRACT 


Article history: 
Received:23 August 2017 
Revised: 27 August 2017 
Accepted: 28 August 2017 


Keywords: 
Deep beam; 
Deformations; 
Deflections; 
Aluminum. 


Many parts of spacecraft, airplanes are made up of 
aluminum, due to its property of less density, which is deep 
in sections. For the analysis of deep parts of any structure, a 
shear deformation theory using the trigonometric sinusoidal 
function in displacement field in terms of thickness 
coordinate is developed to obtain the shear deformation 
effects. The shear stresses are obtained from the use of 
constitutive equations with outstanding accuracy, satisfying 
the zero shear stress at the both, top and bottom of beams. 
Also, the theory not requires the shear correction factor. By 
using the principle of virtual work, the governing differential 
equations and boundary conditions are obtained. A deep 
aluminum beam is assumed subjected to a cosine load for the 
analytical study to show the accuracy of the theory. The 
results are compared with other theories. 


1. Introduction 


Euler-Bernoulli theory disregards the property of the shear deformation and heavy stress 


concentration hence it is applicable for thin beams only. As the theory neglect the shear 
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deformation, and it is more pronounced in the case of deep beams. Hence it is essential to 


address the shear deformation effects. 


Timoshenko [1] considered a prismatic bar to investigate the effect of transverse vibration. This 
theory is known as first-order shear deformation theory (FSDT) and Timoshenko theory in the 
literature. The shear strain is assumed to be constant through the thickness of the beam and thus 


requires shear correction factor. 


Cowper [2] has denied an expression for the shear correction factor for various cross-sections of 
the beam. The discrepancies in the elementary theory of beam and first-order shear deformation 
theories; various higher order shear deformation theories were developed and available in the 
literature for static and dynamic analysis of beams. Krishna Murthy [3], Baluch et al. [4], A 
parabolic shear deformation theory is developed by Bhimaraddi and Chandrashekhara [5] 
assuming a variety of axial displacement. These theories satisfy shear stress conditions on both 


faces of the beam and thus no need of shear correction factor. 


Dahake and Ghugal [6] studied flexural analysis of thick beam having simple supports using 
trigonometric shear deformation theory. Ghugal and Dahake [7,8] given the flexural solution for 
the beam subjected to parabolic loading. Sawant and Dahake [9] developed the new hyperbolic 
shear deformation theory. Chavan and Dahake [10,11] presented clamped-clamped beam using 
hyperbolic shear deformation theory. The displacement and stresses for the thick beam given by 
Nimbalkar and Dahake [12]. 


Jadhav and Dahake [13] presented bending analysis of deep cantilever beam using steel as a 
material. Manal et al. [14] investigated the deep fixed beams using new displacement field. Patil 
and Dahake [15] carried out finite element analysis using 2D plane stress elements for the thick 
beam. Dahake et al. [16] studied the flexural analysis of thick fixed beam subjected to cosine 
load. Tupe et al. [17] compared various displacement fields for static analysis of thick isotropic 


beams. 


In available research, most of the researchers have used steel as a beam material. As many parts 
of the spacecraft, airplane structures are made up of aluminum due to its low weight density. In 
this research, an attempt has been made to analyze the aluminum deep cantilever beam subjected 


to cosine load. 


2. Development of theory 


The assumed beam is made up of homogeneous and isotropic material which occupies in 


O—x-—y-z the Cartesian coordinate system the region: 
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O<x<L; Os y<b; -h/2<z<h/2 
where x, y, z = Cartesian coordinates, 
L, b = length, and width of the beam in the x and y directions, respectively, and 


h = thickness of the beam in the z-direction. 


2.1 Displacement field used 
The displacement field used as follows: 


dw 


uUCx,z)= -2( ae 


) + = [sin — = 60 | and w(x, z) = w(x) (1) 
a h 


Where, uv and w = axial and transverse displacements respectively. 


Plane strain 


2 
h 

ca aL AS La (2) 

ox dx Ki h dx 
Shear strain 

Ou dw NZ 

= —+— =cos— 3 
aa tae - Y (3) 

Stress-Strain Relationships 
o, =Eé (4) 
i = Gy. 


2.2. Boundary conditions and governing equations 

The variationally consistent governing differential equations and boundary conditions for the 
beam under consideration are mentioned below using the expressions for strains and stresses (2) 
through (4) and using the principle of virtual work; it leads to: 


bf i ee (o,62, + T, OY )Ax dz [oa Oowdx =0 (5) 


Where, 6 = variational operator. 


Employing Green’s theorem in Eqn. (4) successively, we obtained the coupled Euler-Lagrange 
equations which are the governing differential equations and associated boundary conditions of 


the beam. 
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4 3 
d‘w _24 , @¢ 


EI = q(x 6 

dx* 5 oa dx? gO (6) 
24 d*w_ 6 dp GA ,_ 7 
tl ag aa a He yO =O (7) 


Moreover, associated consistent natural boundary conditions obtained are as follows: 


Atx=Oandx =L 


2 
Ag 9) cL (8) 


EI—, = 
dx dx 


2.3. General solution of equilibrium equations 


Using Eqns. (6) and (7), we get w(x) and ¢ (x) 


d'w _24d°9 Q(x) 


= 9 
dx =m dx’ EI 0) 


where Q(x) = generalized shear force for beam. 


The second governing Eqn. (7) is also rearranged as follows: 


d’w ado Bo 
dx’ 4 dx’ (10) 


Using Eqns. (11) and (12), we obtained ¢ as: 


dp 45, Ox) 
i a = 11 
dx? eo aE! (1) 


where constants a, and /4 in Egns. (10) and (11) are as follows 


The general solution of Eqn. (11) is obtained as follows: 


= : 12 
P(x) = C, cosh Ax + C, sinh Ax : (12) 


The equation of w(x) is obtained by substituting the expression of ¢(x) in Eqn. (12) and then 


integrating it thrice with respect to x. 
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The solution for w(x) is obtained as follows: 


Ce iCk El 3s a 
El w(x) =| [[[adedededx + - [Ze p\a (C,sinhAx + C,coshax)+C, + C,x+C, (13) 


where, C,,C,,C,,C,,C, and C,are arbitrary constants and can be obtained by imposing natural 


(static) and / or geometric or kinematical end conditions of the beam. 


3. Illustrative example 


To establish the effectiveness of the theory, a numerical example is assumed. For the static a 
rectangular deep cross-section, having span ‘L’, width ‘b’ and thickness ‘h’ of elastic material is 


considered. The following material properties for beam are used. 


Table 1 
Properties of aluminum 6061-T6, 6061-T651 [18]. 


Physical Properties Quantity 


Density 2700 kg/m? 


Ultimate Tensile Strength | 310 MPa 


Modulus of Elasticity 68.9 GPa 


Notched Tensile Strength | 324 MPa 


Ultimate Bearing Strength | 607 MPa 


Bearing Yield Strength 386 MPa 


Poisson’s Ratio 0.33 


3.1. Example: A clamped beam subjected to cosine load 


The beam has its origin at x = 0 and free at x = L. The load acting as shown in Fig. 2, on surface z 


= +h/2 acting in the downward z-direction. 
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TX 
qo q(x) = q og 


Fig. 2. Beam subjected to cosine load. 


Associated boundary conditions for the above beam are: 


2 3, 2 
Apecoenany wy!" a ary 0g) 
dx dx dx dx 


and at fixed end: - = d=w=0atx=0 
X 


The expressions obtained in general form for w(x) and ¢(.x) are as: 


960 °( 7X 1x? | 4920{1x° 1x 
a | OO ae eee ee Ie ee 
(i 1 2L 2L 1 6L 4L 


4 2 2 
moe Gol. a EI °{ eo Xx 7 Ta 


S 
2L 


~1200E1| «GL | a oe bs 
: £2 (-2-j5 +See) 
2 2 
Cie) toi AL at 
be ae 
(x)= fom | 4 —-= gin 7” + sinh Ax—cosh Ax (15) 
BEI. 41 2L 


The stresses and axial displacement are obtained using the above solutions are as follows 
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4. Results 


4.1. Numerical Results 


In the present research, the results for displacements, (axial, transverse and inplane) also, stresses 


are investigated. The results are presented in the non-dimensional form. 


Ebu _ 10Ebh>w _~ bo _ bt 
w=—__, 0,=—, 1T,=—~ 


qh Gi do Ao 


u= 


The shear stresses (7,.) are obtained directly by using constitutive relation and equilibrium 
equation of two-dimensional elasticity and are denoted by (7S*) and (72") respectively. It 
satisfies the stress-free conditions on the top (z=—h/2)and bottom (z=+h/2) faces of the 


beam. 


78 P. Kapdis et al./ Journal of Soft Computing in Civil Engineering 2-1 (2018) 71-84 


Table 2 
Axial Displacement (i ) at (x = L, z = h/2), Deflection (iW ) at (x = L, z =0.0) Bending Stress (0, ) at (x 
= 0, z = h/2), Shear Stresses 75% (x=0.01L, z= 0) and 7%” (x, z= 0) of the Beam for Aspect Ratio 4 
and 10. 


Source Aspect ratio Model a w os Ga ‘a 
Present TSDT -67.5989 6.1819 36.7529 1.8181  -2.7877 
Sawant and Dahake [9] HPSDT -70.024 6.1928 39.8104 2.1609 -4.5581 
Krishna Murty[3] 4 HSDT -71.2291 6.1860 37.2887 1.9004 -2.8916 
Timoshenko [1] FSDT 23.1543 6.5444 22.2081 0.3076 3.7597 
Bernoulli-Euler ETB 23.1543 5.7541 22.2081 — 3.7597 
Present TSDT = -1055.4548 95.8244 176.7877 = 7.7501 = 3.2052 
Sawant and Dahake [9] HPSDT -1061.5175 5.8256 178.4137 8.3023 3.8042 
Krishna Murty[3] 10 HSDT — -1064.5122 5.8255 172.1017 7.8208 3.7523 
Timoshenko [1] FSDT 361.7856 5.8805 138.8010 4.8073 9.3993 
Bernoulli-Euler ETB 361.7856 5.7541 138.8010 _ 9.3993 


Present TSDT 
HPSDT 
HSDT 


FSDT 0.2 


eee 


ETB 


-75 -60 -45 -30 -1 15 30 45 60 75 


u 


Fig. 2. Variation of « through the thickness of the beam for aspect ratio 4. 
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Fig. 3. Variation of « through the thickness of the beam for aspect ratio 10. 
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Fig. 4. Variation of a maximum w of the beam with aspect ratio S. 
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Present TSDT 
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Fig. 5. Bending stress (©, ) through the thickness at (x= 0, z) for aspect ratio 4. 
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Fig. 6. Bending stress (© , ) through the thickness at (x= 0, z) for aspect ratio 10. 
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Fig. 7. Shear stress (7 ,,.) through the thickness at (x = 0.01L, z) obtain using CR for aspect ratio 4. 
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Fig. 8. Shear stress ( 7 ,.) through the thickness at (x = 0.01Z, z) obtain using CR for aspect ratio 10. 
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Fig. 9. Shear stress ( 7 _.) through the thickness at (x = 0.01Z, z) obtain using EE for aspect ratio 4. 
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Fig. 10. Shear stress ( 7 _,.) through the thickness of beam at (x = 0, z) obtain using EE for aspect ratio 10. 


5. Conclusions 


The theory gives reasonable results of this component in the displacement field in proportionate 
with various shear deformation theories. The transverse deflection given is in good agreement 
with that of other higher order shear deformation theories. The flexural stress and its distribution 


across the thickness given by theory are in very good agreement with that of higher order shear 
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deformation theories. Also, the for a cantilever beam with cosine load, transverse shear stress and 


its distribution through the thickness of beam obtained using constitutive relation is in close 


agreement with that of other higher order refined theories; however, it cannot predict the effect of 


stress concentration at the built-in end of the beam. By using the equilibrium equation of two- 


dimensional elasticity, the effect of stress concentration is exactly predicted. 
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